On the equivalence of the freely cooling granular gas to the sticky gas 
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A freely cooling granular gas with velocity dependent restitution coefficient is studied in one di- 
mension. The restitution coefficient becomes near elastic when the relative velocity of the colliding 
particles is less than a velocity scale 5. Different statistical quantities namely density distribution, 
occupied and empty cluster length distributions, and spatial density and velocity correlation func- 
tions, are obtained using event driven molecular dynamic simulations. We compare these with the 
corresponding quantities of the sticky gas (inelastic gas with zero coefficient of restitution). We find 
that in the inhomogeneous cooling regime, for times smaller than a crossover time ti where ti ~ 
the behaviour of the granular gas is equivalent to that of the sticky gas. When (5 — > 0, then ti oo 
and hence, the results support an earlier claim that the freely cooling inelastic gas is described by 
the inviscid Burgers equation. For a real granular gas with finite 5, the existence of the time scale 
t-i shows that, for large times, the granular gas is not described by the inviscid Burgers equation. 

PACS numbers: 45.70.Mg, 47.70.Nd, 05.70.Ln 
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I. INTRODUCTION 

Granular systems, ubiquitous in nature, exhibit a 
wide variety of physical phenomena P, 0. The non- 
equihbrium nature of these systems, characterized by dis- 
sipative particle collisions and external driving, has at- 
tracted a lot of attention from the theoretical point of 
view ^3]. However, analytical progress remains slow. It 
is therefore important to understand completely simple 
models that capture some essential feature of these sys- 
tems. 

A simple, well studied model that captures dissipa- 
tion effects, is the freely cooling gas wherein particles 
at high temperature come to rest through inelastic coUi- 

sion£iJi,Ja M i, 0, [H [H El, Q m, 13 113, d, 

[l9l.l20ll2ir[23.l23|. At short times, the particles undergo 
homogeneous cooling, described by the Half's law |24| . 
wherein the kinetic energy of the system decreases with 
time t as . This behaviour has recently been veri- 
fied experimentally [25j . Beyond a crossover time tc, the 
system crosses over to the inhomogeneous cooling state, 
characterized by clustering, where the energy decays ac- 
cording to a different power law ~ t""^. The exponent ip 
is dependent on dimension. For instance, i/i = 2/3 in one 
dimension and i/> w 1 in two dimensions. 

In one dimension, it has been possible to make con- 
siderable analytic and numerical progress. In the spe- 
cial case when the coefficient of restitution (henceforth 
denoted by r) is zero, i.e. particles stick on collision, 
the problem is exactly solvable. The problem can be 
mapped on to the description of shocks in velocity field 
in the inviscid Burgers equation 4]. Burgers equation 



being exactly solvable in one dimension l2q|, the exact 
expressions for different quantities like velocity distribu- 
tion, mass distribution, energy decay as a function of 
time are known. This limit (r = 0) will be referred to 
as the sticky gas. In this paper, we will refer to the case 
< r < 1 as the granular gas. 

Extensive simulations of the granular gas in one di- 
mension showed that, beyond the crossover time tc, the 
energy decay in the granular gas is identical to that of 
the sticky gas Not only is the exponent if) the same, 
all the curves for different r's collapse on top of each 
other. In addition, the velocity distribution was seen to 
be independent of r. Finally, the coarsening length scale 
C{t) scaled with t identically for both the granular and 
the sticky gas. It was thus concluded that the long time 
behaviour {t > tc) of the granular gas is described by 
the Burgers equation. This analogy was generalized to 
higher dimensions in Ref. [6j] , though its applicability re- 



mains unclear [17 1 
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To make this equivalence in one dimension stronger, 
the two point correlations of the granular and the sticky 
gases need to be compared. In a recent paper we 
showed that, when r depends on the relative velocity, 
there is a new crossover time ti (with ti > tc), beyond 
which the coarsening properties of the granular gas differ 
from that of the sticky gas or the Burgers equation. This 
was shown numerically by measuring the spatial density- 
density and velocity-velocity correlation functions. In 
this paper, we focus on the intermediate time tc < t < ti, 
and show that the correlation functions of the granular 
gas compare well with that of the sticky gas, obeying 
what is known as the Porod law. In Sec. [Tll we review 
this law in the context of coarsening systems, and point 
out some non-equilibrium systems where it is known to 
be violated. 

Our study involves event driven molecular dynamics 
simulations of the granular and the sticky gases on a ring. 
The precise definition of the model and the definition of 
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the statistical quantities that we measure are presented 
in Sec. IIIII The quantities that will be studied are den- 
sity distributions, density-density and velocity-velocity 
correlation functions, and size distribution of empty and 
occupied clusters. 

Modeling the coefficient of restitution r is crucial in 
defining the model and understanding cooling granular 
gases. Experimentally, it is known that r tends to 1 (elas- 
tic limit Vwhen the relative velocity of collision decreases 
to zero psi . [29l [30l |. while for large relative velocities r 
tends to a constant. This introduces a velocity scale 5 in 
the problem, which consequently gives rise to a new time 
scale ti. As in Ref. [131, we model r as 

r = (l-ro)exp(-|i;,ei/'^r)-t-ro, (1) 

where 5 represents a velocity scale, below which colli- 
sions become near elastic, i.e, when the relative velocity 
Vrei ^ 5, then r — > 1 (elastic collision). For Vrei ^ 5 
the collisions are inelastic with r — > tq. The parame- 
ter tj decides how abruptly this crossover occurs. For 
example, for cr ^ oo, r versus Vrei is a step function 
while for finite cr it is a smoother curve. For cr < 1, the 
curve has a cusp as it approaches r — 1. The exper- 
imental curves [28| suggest that the power a can vary 
over wide range. Within the framework of viscoelastic 
theory, a power of cr = 1/5 has been predicted [3l|, . 
Apart from that, studies of velocity distribution in freely 
cooling granular gases with velocity dependent restitu- 
tion coefficient have been done using techniques different 
from that in this paper, namely direct simulation monte 
carlo and viscoelastic theories [sl, [11] . With a ^ oo and 
5 held finite, our problem becomes same as the granular 
gas problem studied earlier in 

The existence of a new crossover time ti was briefly 
shown in our earlier publication [23]. In Sec. IIV( we 
present simulation results, showing explicitly the exis- 
tence of ti for different 5 and a and determine its de- 
pendence on 6 and a. We then proceed to do a detailed 
comparison of the early time behaviour of the granular 
gas and the sticky gas in Sec. [V] We do the comparison 
on the basis of data for the various statistical quantities 
characterizing coarsening in these systems. In Sec. IVI) 
the two point correlations are related to the empty and 
occupied cluster distributions through a mean field ap- 
proximation. This derivation leads us analytically to the 
Porod law for these coarsening systems. Conclusions and 
discussions are presented in Sec. IVIII 



II. THE POROD LAW 

The freely cooling granular gas is an example of what 
is known as a coarsening or a phase ordering system. It 
is well known that in systems freely relaxing to ordering 
states, there exists a length scale C{t) that increases with 
time [35| . In addition, for usual phase ordering systems, 
the presence of a dominant C{t) results in a robust scaling 



law called the Porod law [s^, [s^ . For scalar order pa- 
rameters, the Porod law states that in c?-dimensions, the 
scaled structure function S{k,t)/C'^ ^ (kC)^^ for large 
kC, with — d+1. Hence, in one dimension 9 — 2. In the 
above, S{k, t) is the Fourier transform of the two point 
correlation function of the order parameter. We shall re- 
fer to systems having the above features as clean phase 
ordering systems. The essential fact is that, in clean 
phase ordering systems the relative fluctuation with re- 
spect to the ordered state is nominal. For example, in the 
case of a ferromagnetic Ising model, a finite system at in- 
finitely large time shows magnetization per site peaked 
at two possible values of spontaneous magnetization. 

A few non-equilibrium systems, e.g., p articles sliding 
on a randomly fluctuating surface |37l. ISa. Is^ . active ne- 
matics namely, agitated granular-rod monolayers or films 
of orientable amoeboid cells [i^l, show a different kind 
of phase ordering; the main characteristic of which is 
undamped relative fluctuations even in the thermody- 
namic limit. Hence, it is called fluctuation dominated 
phase ordering. In such systems strong fluctuations sus- 
tain in time without losing macroscopic order. The or- 
der parameters have broad probability distributions even 
in the thermodynamic limit indicating strong variation 
of the order parameter in time. Unlike usual phase or- 
dering, there need not be sharp interfaces distinguishing 
one phase from the other. The ensemble-averaged spa- 
tial correlation function shows a scaling form in \r/£{t)\, 
but unlike clean phase ordering systems, it exhibits ei- 
ther a cusp or a power law divergence at small values 
of \r/£{t)\. Hence the scaled structure factor varies as 
S{k,t)/C'^ ~ (kC)-" with 9 < 2. Thus, measuring 61 is a 
reliable test to identify deviations from the Porod law. 

We recently showed that the freely cooling granular 
gas at large times {t > ti) also shows fluctuation domi- 
nated phase ordering [23|. The spatial correlation func- 
tion shows a power law divergence for small values of 
\r/C{t)\. This in turn indicates violation of the Porod 
law in the density and velocity structure functions with 
9 « 0.8 [23|. The Porod law violation here indicates the 
presence of power law distributed clusters. In other words 
the phase ordering is not clean. Earlier, in a slightly dif- 
ferent model in one dimension with velocity dependent 
restitution coefficient [4lj, it was found that the coars- 
ening clusters are not long lived and they start breaking 
up at late times. That result seems to have similarity 
to what we found in at late times t > ti. But in- 
terestingly we find that at late times {t > ti) there are 
two sub- regimes namely t2 > t > ti and t > t2- Here 
t2 ~ S^^ ^ is the very large scale beyond which all 
collisions become elastic and the system again becomes 
homogeneous. On the other hand, the regime t2 > t > ti 
is the new and nontrivial regime of "unclean" ordering 
that we have found [l^l and are highlighting here. 

For Burgers equation, it was shown numerically in 
Ref. [l^l [also see Sec.|V], that the ordering is clean phase 
ordering, such that Porod law is obeyed. In this paper 
we examine the structure factor for the granular gas for 
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intermediate times and show that Porod law is obeyed, 
Ukc in the sticky gas. 

III. THE MODEL 

In this section, we define the model and the different 
physical quantities that are measured. Consider TV point 
particles of equal mass on a ring of length L. Initially, 
the particles are distributed randomly in space with their 
velocities drawn from a normal distribution. Choose the 
reference frame in which the center of mass momentum 
is zero. The particles undergo inelastic, momentum con- 
serving collisions such that when two particles i and j 
with velocities Ui and Uj collide, the final velocities u'^ 
and u'j are given by: 

where the relative velocity dependent restitution coeffi- 
cient r is given by Eq. ([1]). For the sticky gas, 6 = and 
rg = 0, such that r = for all collision velocities, and 
the particles coalesce on collision. 

For both the systems, the granular gas (0 < tq < 1) 
and the sticky gas (r = 0), we have done event driven 
molecular dynamics simulations with system size L = 
20000 in units of the mean inter particle spacing at time 
t = 0. The particle density is set to 1 throughout. Peri- 
odic boundary conditions were imposed. For the granular 
gas, simulations were performed for ro = 0.1, 0.5 and 0.8 
in the time regime in which C{t) <C L. Thus, formally the 
sequence of limits are L — > oo first and then t oo. In 
such a regime, the system continually coarsens, and hence 
we refer to it as a "phase ordering system" . In this time 
regime, simulations were done for various different values 
of a (namely 3,4,5,10 and oo), S (namely 0.001, 0.002 
and 0.004). We found no qualitative variation of our 
main results with varying values of ro,cr and S. There- 
fore, unless otherwise specified, the data that is presented 
for the granular gas corresponds to ro = 0.5, a = 3.0, and 
S = 0.001. 

It is well known that for (5 = 0, event driven molecu- 
lar dynamics simulations cannot go beyond the inelastic 
collapse. A similar difficulty is faced when 5 is finite, but 
CT 0. For example, doing simulation with realistic a 
like 1/5 (suggested from viscoelastic theories) or 1/4 ( for 
plastic collisions at large velocities), is out of our reach. 
At such small a the shape of the r versus Vrei curve has 
a cusp as Vrei — > 0. The particles practically do not see 
an elastic regime and the simulation run into difficulties 
like particle overlap and inelastic collapse. Nevertheless, 
from our studies of (p^) versus t (fig. 3), we see that with 
decreasing a the "clean phase ordering regime" gets ex- 
tended further in time. So we hope this trend would 
continue even for smaller a that our simulation could 
not access. Also, the exponents are independent of the a 
values that we have studied, and we hope this trend will 
hold for smaller a too. 



The primary difference between a sticky gas and an 
granular gas is that in the former particles coalesce to 
form larger and larger mass clusters, while in the latter 
masses never coalesce. As a result, "clustering" phenom- 
ena in granular gas can be meaningfully quantified only 
if we define a local coarse-grained mass density p. To do 
so, first divide the ring into equally sized N boxes. For 
the granular gas, pi is defined to be the total number of 
particles in the ith (i = 1,2,3,..., N) box. The coarse 
grained velocity Vi is defined as the mean of the velocities 
of all the particles in box i. 

We now define the statistical quantities that are mea- 
sured for the granular gas. The extent of clustering can 
be measured in terms of the variation of the density dis- 
tribution P{p, t). It is the probability of finding a spatial 
box with density p at time t. The first moment {p){t) 
of the distribution P(/9, t) is constant at all times as the 
total mass is a conserved quantity. The second moment 
{p^) (t), is defined in our simulation as {J^i pi^) /N, where 
the sum is over all the boxes at time t, and the average 
( ) is over initial conditions. 

Other measures of clustering are the cluster lengths 
and gap lengths. The distribution 0{x,t) is the number 
of clusters with exactly x successive occupied boxes at 
time t, normalized by the total number of occupied clus- 
ters. Likewise, the distribution E{x,t) is the number of 
clusters with exactly x successive empty boxes at time 
normalized by the total number of empty clusters. 

The other quantities of interest in this paper are spa- 
tial correlation functions. The equal time, two-point 
density-density spatial correlation function is defined as 
Cpp{x,t) = {pi(t)pi+^(t)), where pi and pi+x are densi- 
ties at box i and box (i + x) respectively. The structure 
function Spp{k,t) is the Fourier transform of Cpp{x,t). 
Similarly, the equal time, two-point velocity-velocity 
spatial correlation function is defined as Cvv{x^t) = 
{vi{t)vi+xit)), where Vi and Vi+x are the coarse grained 
velocities of box i and box (i -I- x) respectively. The 
structure function Svy{k,t) is the Fourier transform of 

In comparison to the granular gas, as particles coalesce 
in a sticky gas, they form bigger and bigger mass clusters 
localized at points in space, such that their total number 
steadily decreases with time. Therefore, without coarse 
graining, one can measure mass and velocity distribution 
functions, and the corresponding two-point correlation 
functions in space. But, to maintain uniformity with 
the granular gas, we will continue to use coarse grained 
densities and velocities for the sticky gas. 



IV. THE NEW CROSSOVER TIME h 

In this section, we show that the introduction of a ve- 
locity scale 6 introduces a new time scale ti, beyond 
which the granular gas deviates from the usual clean 
phase ordering as seen in the sticky gas. We first show 
the different coarsening regimes by means of a space-time 
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FIG. 1: The spatial density distribution is shown for different 
time regimes: (A) t = - 7 (t < tc) {B)t = 600 - 607 
{tc <t <U) and (C) t = 8000 - 8007 (t > ti). The darker 
regions correspond to higher densities. The data is for the 
granular gas with r = 0.5, 5 = 0.004 and a — oo. For these 
values of the parameters, ti ~ 1500. 



density plot. In Fig. [1] the density distribution is shown 
for three different time regimes with the darker regions 
corresponding to higher density. 

In Fig[IJA), we show a snapshot of the density profile 
of the system in a finite portion of space for few succes- 
sive time steps, for very early times, t ^ tc- We see that 
the system is fairly homogeneous. In Fig[TjB), a simi- 
lar plot of the density profile is shown for intermediate 
times tc < t < ti. The picture looks like any ordinary 
clean phase ordering system with steady growth of iso- 
lated large density lumps (deep dark shaded), along with 
large empty gaps (white) separating them. If the clean 
phase ordering continues, then at later times one would 
expect one or two very huge lumps remaining. Instead, in 
Fig.[ljC), we find grey shaded spatially extended patches 
reappearing, along with large white gaps, implying frag- 
mentation of some large clusters and growth of smaller 
clusters. This is the signature of fluctuation dominated 
phase ordering. 

The above picture can be quantified by measuring the 
second moment of the density distribution (p^). For the 
sticky gas, it is known that (p^) ~ i^/^ for t 1. Does 
(p^) in the granular gas behave as that of the sticky 
gas? Since there are three parameters, cr, d and rg in 
the model, we need to carefully address the case with 
respect to varying these parameters. 

Figure [2] shows the variation of {p'^){t) with time t for 
varying 5, with a = 3.0 and rp = 0.5 fixed. The trends 
are the same for other values of a. For early times tc < 
t < ii, (p^) ~ t^/^ and completely overlaps with the same 
data for the sticky gas (shown by a line). For late times 
t^ti, {p^) deviates from the sticky gas behaviour, the 
crossover occurring later for smaller 5 (see Fig. [2])- The 
quantitative dependence of the crossover time ti on 5 can 
be obtained by collapsing the different curves by scaling. 
The collapse is excellent (see inset of Fig. [5]) when the 




FIG. 2: The temporal variation of the second moment of den- 
sity, (p^) for the granular gas, is shown for different values of 
5 with a = 3 for all the cases. The data is for the granular gas 
except for the dotted line which corresponds to the sticky gas. 
The dotted line varies as t^^^ . The data for the granular gas 
coincides with that of the sticky gas for intermediate times, 
before deviating. The inset shows data collapse for various 5 
when the variables are scaled as in Eq. ([3}. 



different curves and t are scaled as: 



(3) 



Hence, we conclude that ti ~ 6^^ . 

The crossover time ti also depends on a. Figure [3] 
shows the variation of (p^) with time t for different values 
of cr, keeping 5 fixed at 0.001 and tq = 0.5. As earlier, 
for intermediate times the behaviour of (p^) mimics that 
of the sticky gas, while for large times, there is deviation 
beyond a scale ti. The crossover time increases with 
decreasing a. When a = oo, ti is finite but nonzero. 
Thus, we conclude that for all values of cr, the crossover 
time ti is finite and nonzero. 

Finally, we have checked that for fixed 5 and a, curves 
of (p^) do not vary with tq, thus implying that ti does 
not depend on tq. 

The coarsening behaviour for times t ^ ti was dis- 
cussed extensively in Ref. pTj . In this paper, we re- 
strict ourselves to the regime tc < t < ti, and compare 
the coarsening behaviour of the granular gas to that of 
the sticky gas. We bring out the similarities of the two 
through extensive simulations, the results of which are 
given below. 



V. COMPARISON OF THE GRANULAR GAS 
TO THE STICKY GAS 

In this section, we present results for the granular gas 
and the sticky gas for the intermediate times tc < t < ti. 
We show that they compare very well, and thus conclude 
strongly in favor of the equivalence of the two in this time 
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FIG. 3: The temporal variation of the second moment of 
density, {p^) for the granular gas, is shown for different values 
of a for 5 — 0.001. The curve with power (p^) ~ t^^^ overlaps 
with all the data for tc < t < t\. The crossover time ti 
increases with decreasing cr, and is finite for a = oo. 

regime. The following quantities are studied: the den- 
sity distribution t), the occupied cluster distribution 
0{x,i), the empty interval distribution E{x,t), density- 
density correlations Cpp{r, t) (in real space) and Spp(k, t) 
(in Fourier space) and the velocity-velocity correlations 
Cm{r,t). For the sticky gas, the scalings of these dis- 
tribution functions are known exactly through the exact 
solution in ^ and earlier studies However, some 

of the results, for example the velocity-velocity correla- 
tions, are complicated expressions whose behaviour is not 
easily visualized. In addition, no plots are provided for 
ready reference. For completeness and future reference, 
we present numerical data for the sticky gas generated 
from event driven simulations. The results are, of course, 
in agreement with d, [ll] . 

A. Density Distribution function 

The mass or density distribution for the sticky gas is 
a power law with a cutoff increasing with time, and an 
amplitude that decreases with time For the sticky 
gas, it has the scaling form 

where the scaling function /i (z) ~ 1/ ^/z when 2^1 and 
/i(z) when z ^ 1. Thus, for densities much smaller 
than the cutoff, P{p,t) ity/p)~^- The simulation re- 
sults for the sticky gas are shown in Fig. |4] for different 
times t. The data shows a time dependent cutoff, and 
power law with exponent —1/2. The scaling behaviour 
is consistent with Eq. (jll (see inset of Fig.[4|). 

Now, consider the granular gas for tc < t < ti. The 
variation of the density distribution P(p, t) with p for 
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FIG. 4: The variation of P{p,t) with p is shown for different 
times for the sticky gas. The straight line has an exponent 
— 1/2. The inset shows the data collapse when scaled as in 
Eq. gl). 
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FIG. 5: The variation of the density distribution P(p, t) with 
density p is shown for different times (in the early time regime) 
for the granular gas. The straight line has an exponent —0.5. 
The inset shows the data collapse when scaled as in Eq. Q. 



different times is shown in Fig. [S] The distribution is a 
power law with the cutoff increasing with time and the 
amplitude decreasing with time. The cutoff scales exactly 
as in the case of sticky gas, i.e., Pmax ~ t^^^- A good data 
collapse (shown in the inset of Fig. [5]) is obtained when 
the different data are scaled as in Eq. ([4]). 

We note that if (5 ^ the behavior of P{p, t) as in 
Fig. [5] would continue indefinitely in time just as in sticky 
gas. But, in a more realistic granular gas model with fi- 
nite (5, for t 3> ii, the behavior of P{p,t) completely 
changes and a new scaling function emerges asymptoti- 
cally as was detailed in Ref. [27i] . 
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FIG. 6: The variation of the empty interval distribution 
E{x,t) with separation x is shown for the sticky gas at dif- 
ferent times. The inset shows the collapse of the curves when 
scaled as in Eq. ([5]). 



B. Empty and Occupied cluster Distribution 
Functions 

A suitable measure of clustering in space is the empty 
gap distribution E{x,t). For the sticky gas it is known 
Hi that 

E{x,t)^^h[^), (5) 

where the scaling function f2iz) — ^ constant when z — > 
and f2{z) — *■ when z ^ 1. Thus, the largest empty clus- 
ter sizes Xmax ^ t^^'^ ■ The simulation data for E{x, t) for 
the sticky gas is plotted in Fig. [5] for various times. The 
increasing cutoff points to the steady increase in the size 
of empty clusters, as is usual in phase ordering systems. 
The decreasing amplitude of the distribution for small 
sizes indicate that small empty gaps are steadily disap- 
pearing from the system with time. In the inset of Fig. [6l 
the data collapse is shown following Eq. ([5]). 

For the granular gas, the scaling and shape of the scal- 
ing function for E{x,t) matches well with that of the 
sticky gas. In Fig.[7|we present the data for various times 
and in the inset the scaling collapse using Eq. ([5]). The 
collapse is excellent, showing the equivalence of the gran- 
ular gas to the sticky gas as far as E{x, t) is concerned. 

It turns out that in the granular gas for t ^ ti, the 
above clean phase ordering behaviour of E{x,t) dramat- 
ically changes and a power law distribution for the scal- 
ing function /2(z) appears with a large negative power 
(w 2.2) [131, signaling a crossover to fluctuation domi- 
nated phase ordering behaviour. 

The occupied cluster distribution 0{x,t), has a very 
different form and evolution in comparison to E(x, t) 
discussed above. In the absence of any analytical guid- 
ance on this quantity from earlier works on sticky gas. 
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FIG. 7: The variation of the empty interval distribution 
E{x,t) with separation x is shown for the granular gas at 
different times. The inset shows the collapse of the curves 
when scaled as in Eq. ([Sjl. 
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FIG. 8: The variation of the occupied cluster length distribu- 
tion 0{x,t) with separation x, is shown at different times for 
the sticky gas on a linear-log plot. Inset: A similar plot for 
the granular gas. The times are t — 200(top curve), 400,800 
and 1600(bottom curve). 

we present directly what we find in our simulation. The 
data for 0{x,t) is shown in Fig. [S] We see that large 
occupied clusters are extremely rare, and within size 3 
the distribution fall to very low value. The distribution 
function 0{x, t) is exponential with x, and its width de- 
creases marginally with time. Thus, we conclude that 
occupied clusters are highly localized in space and do 
not vary much in time. 

In the inset of Fig. [SI we show 0{x, t) for the granular 
gas in the early time regime. The curves are strikingly 
similar to the sticky gas. There is an exponential decay 
over very short sizes, and clusters beyond size 9 are very 
rare. Again there is a marginal decrease in width of the 
distribution with time. Just as the sticky gas, this implies 
extremely localized occupied clusters, whose width do 
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FIG. 9: The density-density correlation function Cpp as a 
function of a; for the sticky gas. Inset: Shows the data collapse 
when scaled as in Eq. ([6}. 



FIG. 10; The density-density correlation function Cpp as a 
function of x for the granular gas. Inset: Shows the data 
collapse when scaled as in Eq. ((Bj. 



not vary appreciably with time. 

We note that, in sharp contrast to the above scenario, 
the occupied clusters do not remain localized at large 
times t ^ ti, for the granular gas. In fact 0{x,t) also 
decays as a power law like E{x, t) with the same negative 
power (w 2.2) as was shown in Ref. [27| . 

C. Density-density correlation functions 

The two point spatial correlation functions give a very 
vivid picture of the type of phase ordering that the sys- 
tem undergoes, and we discuss these in this and the next 
subsection. The data for two point density-density cor- 
relation function Cpp (defined in Sec. IIIip of the sticky 
gas is plotted in Fig. [31 The correlation function drops 
from a very large value at a: = to a value below 1 and 
then asymptotically increases to 1 as a; 00. The ini- 
tial drop is because near every particle there is a density 
depletion zone created by coalescence of clusters. The 
asymptotic approach is expected as Cpp (p) for large 
X and (p) = 1. For a phase ordering system [s^ the 
scaling form 

Cppix,t)^f3{x/C), (6) 

is expected, and we show the data collapse in the inset of 
Fig. [HI For small scaled distance x/£, the scaling function 
/a increases linearly, implying Porod law [SG^ to be valid. 

The corresponding data for Cpp for the granular gas at 
intermediate times is shown in Fig. 1101 We find a striking 
similarity between Fig. [TO] and Fig. [5| (for the sticky gas). 
The scaling function shown in the inset of Fig. [TOj is very 
similar to the inset in Fig. [9l 

To find the structure function Spp{k,t) which scales 
in the conventional fashion one has to treat the data for 
Cpp{x,t) carefully. Firstly, the value at a;/£ = for the 



scaled Cpp{x, t) do not form a part of the scaling function 
/a and have to be carefully adjusted by hand so as not to 
introduce undesirable distortions in fc-space. This proce- 
dure has been discussed in details in [33| , and we do the 
same in this case. Secondly, 73(2:) does not decay to zero 
at z ^ 00 and well behaved Fourier transform cannot 
be found. So for this special case of the density-density 
correlation function for both the sticky and the granular 
gases, we define Spp{k,t) to be the Fourier transform of 
1 — Cpp{x, t). The latter real space function nicely decays 
to zero at large x and thus Spp{k, t) is expected to be well 
behaved. The data thus obtained are shown in FigfTTjfor 
the sticky gas, and in Fig. [12] for the granular gas (for 
small times). They compare very well, and in the insets 
of Figs. [TT] and [TO] the scaling collapse of the data are 
shown according to 

Spp{k,t)=Cg:i{kC{t)). (7) 

For large kC we see that for both the sticky and the 
granular gas the scaHng function 53(2;) ~ for z ^ 
1, i.e obey Porod law, as is expected for a clean phase 
ordering system. 

For times t ^ <i there is a shift to non-Porod be- 
haviour for the granular gas as discussed in Ref. [13]. 
The functional form of 53(2:) completely changes in that 
regime. 



D. Velocity- velocity correlation functions 

The velocity-velocity correlation function Cm is ex- 
pected [131 to scale as ~ {vt/ CY fiix/ C), where vt is the 
typical velocity which scales as and the factor 1/L 

accounts for probability of a box being occupied. The 
scaling function fi{x) decays to zero for large x. Putting 
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FIG. 11: The variation of Spp with k for the sticky gas is 
shown for different times. The inset shows the data coUapse 
when scaled as in Eq. ([7]). The straight hne has an exponent 
-2.0. 



FIG. 13: The figure shows the variation of Cvv{x,i) with 
X for different times. The data is for the sticky gas. Inset: 
shows the data collapse when scaled as in Eq. 




,1 



k 



FIG. 12: The variation of Spp with k is shown for the granular 
gas at early times. The inset shows the data collapse when 
scaled as in Eq. ([7}. The straight line has an exponent —2.0. 
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X 

FIG. 14: The figure shows the variation of Cw (x, t) with x for 
different times for the granular gas. Inset: shows the curves 
when scaled as in Eq. (|8]). The higher curves correspond to 
smaller times. At large times, the data approaches the solid 
line, an exponential obtained by fitting the scaling region. 



in the t dependence for vt and £, we obtain 

C™(x,t) - ^/4(a;//:). (8) 

In Fig. [T31 the data for Cvv{x,t) at different times is 
shown for the sticky gas. The inset shows the data col- 
lapse when scaled as in Eq. ([5]). The scaling function 
fi{z) decays linearly for small z and goes to zero for 
large z. The linearity is again a clear signature of the 
Porod law. 

The corresponding data for Cm{x,t) for the granular 
gas is shown in Fig. [T31 There are corrections to scaling 
when compared to the corresponding data for the sticky 
gas (see inset of Fig. [T4)) . A possible reason is that the 
times for which the measurements are taken for the gran- 
ular gas are a factor of ten smaller that of the sticky gas. 
However, it is seen (see inset of Fig. fH)) that, for increas- 



ing time, the x/ C at which deviation from scaling occurs 
decreases. The scaling function tends towards the curve 
3.0exp(— 1.25a;/£) (drawn as a solid line in the inset of 
Fig. [13])). Thus, again a linear decay is obtained near the 
origin, signifying Porod law. The scaling functions in the 
inset of Figs [13] and [M] differ by a factor of 6 for small 
x/ C This is probably a consequence of the coarsening 
rule that we have used. 

We also note that the structure function Svv{k,t), ob- 
tained as the Fourier transform of Cyy{x, t), should scale 
as 

S.,yik,t) ^ ^g^ikC). (9) 

However, due to corrections to scaling for small x/C, we 
cannot obtain good data for large fc, and hence we do not 
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FIG. 15: A schematic diagram showing configurations that 
contribute to pk when fc = 2n — 1 is odd (A) and when k = 2n 
is even (B). The shaded regions represent boxes with nonzero 
density p. The even k correspond to the case when both 
the sites and x are occupied and contribute to the density- 
density correlation function Cpp{x,t). 

show the structure functions as plots. 

VI. THE INDEPENDENT INTERVAL 
APPROXIMATION 

In this section, starting from the empty and occu- 
pied interval distributions, we derive the Porod law for 
density-density correlations via a mean field approxima- 
tion, often referred to as independent interval approxi- 
mation [i^l ■ This approximation was used in Ref . [27'| to 
connect the interval distributions to the density correla- 
tions for times i ^ fi. Here, we first give the derivation 
of the result used in Ref. [23| , and then apply it to times 
tc <t <ti. 

The equal time density-density correlation function 
Cpp{x) — {pi{t)pi^x{i)) is the expectation value of the 



product of densities of boxes separated by spatial dis- 
tance X. This quantity can be estimated using an ap- 
proximation which is discussed below. Let pk{x) denote 
the probability that given that site is occupied, there 
are exactly 1 -I- [k/2] occupied intervals and [{k -I- l)/2] 
empty intervals, where [. . .] denotes the integer part. An 
example of configurations contributing to p„ is shown in 
Fig. fTSTA) for odd fc = 2n — 1 and Fig. [TSTB) for even 
k ~ 2n. Thus, po{x) will correspond to the situation 
when is occupied and there is continuous set of occu- 
pied sites up to a;, while pi{x) would correspond to a set 
of occupied sites starting from 0, followed by an empty 
interval extending beyond x. Since Cpp{x) is approxi- 
mately equal to density square times the probability that 
there is a nonzero density at a; given there is a nonzero 
density at 0, we obtain (setting p = 1), 

oo 

P2n{x), (10) 

ri=0 

where P2n{x) is the probability of having exactly n empty 
gaps between and x with and x being occupied. 

The probabilities Pn{x) may be expressed in terms of 
the gap distributions E{x,t) and 0{x,t). Let 

/"OO 

Q{x,t) = / dx'0{x',t), (11) 

J X 

F{x,t) = / dx'E{x',t). (12) 

J X 

where Q(x,t) is the probability of getting an occupied 
cluster of size greater than x. Similarly F{x,t) is the 
probability of getting an empty cluster of size greater 
than x. 

Approximating joint distributions p2n(a^) s^nd p2n-i{x) 
by products of individual distributions of the intervals 
(namely E, O, Q and F), one can easily verify (following 
Fig. [15] as guideline for the odd and the even cases), that 



2n-3 



P2n-l{x) = a dxi...dx2n-lQ{xi)F{x-X2n-l) Yl i^i^j+i ^ '^3)C{Xj+2 - Xj + i)] , 71=1,2,... (13) 

j = l,3,... 

/2n-3 
dxi . . .dx2nQ{xi)E{x2n - X2n-l)Q{x - X2n) Y\. \-^(^ 1 + ^ ~ ^J^^ix j+2 ~ Xj+i)] , n=l,2, (14) 

j = l,3,... 

Po{x) = a dx'Q{x'), (15) 



I 

where a is a normalization constant that will be deter- we obtain 
mined shortly. 

Equations P^ - P^ simplify considerably if we take _ - - ^ ~vn_i 

Laplace transforms. Denoting the Laplace transform of P2n-i(s) = aQF[EO) , n=l,2,... (16) 

a function f{x) by /(s), where f{s) = f{x)e-'''dx, p2n{s) = aQ^E{Edy'-\ n= 1,2,... (17) 
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Po{s) 



(18) 



where (a;)o = dxxO{x), is the mean length of an oc- 
cupied cluster. The distribution Q{s) can be expressed 
in terms of the probability density 0{s) as 



Qis) 



1 - 0{.s) 



(19) 



and likewise for F (with O replaced in Eq. [12] by 
E). The constant a can now be determined from 
the condition X]^o^"('*) ^ "^tiich follows from 

^^qP„(.t) = 1. We immediately obtain a — l/{x)o- 
Equations HU]) , (US]) , (HH]) and ^ give 



pp(s) = - - 



1 [l-0(.)][l-^(.)] 



s^[l - E{s)0{s)] 



(20) 



We now need the form of 0{s, t) and E{s, t). The occu- 
pied cluster size distribution 0{x,t) has no dependence 
on the length scale C{t) (see Fig. |8]). The distribution is 
very close to a pure exponential and we take it to be of the 
form 0{x,t) « exp(— a;/a) where a is independent of 
time. Then 0{s,t) « (as + 1)^^. The empty interval dis- 
tribution E{x, t) is a constant for small x and goes to zero 
for X ^ Lt (see Figs. [Bland [T]). We approximate it by a 
step function: E{x, t) ~ for x < Ct, and E{x, t) ~ 
for X > Cf Then, E{s,t) « (s/:t)-i[l - exp(-s£f)]. 

We are interested in the limit s — > 0, £ — > oo keeping 
sC a constant. For verifying Porod law, we are further 
interested in the limit sC ^ 1. Substituting the assumed 
forms for E and O in Eq. I|20p and expanding for large 
s£, we obtain 



Co 



1 



c 



sC > 1. 



(21) 



We thus obtain the Porod law behaviour by using the 
cluster distributions. The same independent interval ap- 
proximation gives consistent results between the cluster 
distributions and the two point correlation functions in 
late time regime t > ti [23| , where there is a violation of 
Porod law. 



VII. SUMMARY AND CONCLUSIONS 

In this paper, we studied the density distribution, 
the empty and occupied cluster distributions, the spa- 
tial density-density correlations and the velocity- velocity 
correlations, in a granular gas undergoing inelastic col- 
lisions in one dimension. The restitution coefficient was 
modeled as being velocity dependent, with the collisions 
becoming nearly elastic when the relative velocity is less 
than a velocity scale 5. The velocity scale 5 introduced 
a new time scale ii into the problem. In this paper, 
the dependence of ti on the different parameters was 
found. Due to the existence of ti, and the crossover 
scale tc between the homogeneous and inhomogeneous 
cooling regimes, there are three different time regimes in 
the problem: (1) t < tc, during which particles undergo 
homogeneous cooling, (2) tc < t < ti, referred to as in- 
termediate times, and (3) t > tc referred to as late times. 
In this paper, the focus was on the intermediate times. 

For intermediate times, we compared the different cor- 
relation functions of the granular gas with that of the 
sticky gas. We found that there is an excellent match 
between the two. Coarsening in both the cases is gov- 
erned by the Porod law. These results thus support the 
claim that the Burgers equation is the correct description 
for the granular gas 5]. 

For late times, this equivalence breaks down [23|. The 
coarsening in the granular gas violates Porod law. The 
occupied cluster distribution and the empty cluster dis- 
tributions differ significantly from intermediate times by 
developing into power laws. The existence of clusters 
of all sizes leads to fluctuation dominated coarsening. 
A continuum equation describing this regime is missing, 
and it would be interesting to find one in future. 
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